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ON A NONLINEAR FLUX LIMITED EQUATION 
ARISING IN THE TRANSPORT OF MORPHOGENS 

F. ANDREU, J. CALVO, J. M. MAZON AND J. SOLER 



Abstract. Motivated by a mathematical model for the transport of morphogenes in biological systems, 
we study existence and uniqueness of entropy solutions for a mixed initial-boundary value problem asso- 
ciated with a nonlinear flux— limited diffusion system. From a mathematical point of view the problem 
behaves more as an hyperbolic system that a parabolic one. 
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(_; ' 1. Introduction 

-)— » 

■ ' The aim of this paper is to analyze the mixed initial-boundary value problem associated with a nonlinear 

flux-limited reaction-diffusion system 

dii 

= {a{u,u,))^- f{t-T,u{t,x))u{t,x)+g{t,u{t,x)), in ]0,r[x]0,L[ 



dt 

-a{u{t, 0),u^{t, 0)) == /3 > and u{t, L) = on t e]0, T[, 

u{0,x)=uo{x) in xG]0, L[, 



> 

o : (1-1) 

in 

[~>^ . being 



where the boundary conditions must be interpreted in a weak sense to be precised, and the functions / 
and g are nonlinear with respect to u and depend on u through a coupled system of ordinary differential 
equations. This problem arises in the modelization of the transport of morphogens and the parameter r 



S . represents a delay in the process of signalling pathway cell internalization. 

The nonlinear diffusion equation 

(1.2) ^ = (aKu.)), 

was introduced in different contexts as an alternative to the linear diffusion equation with the ideas of 
limiting the flux and reproducing a system with finite speed of propagation. The flux-limited type equations 
were motivated previously in [3T], but they were firstly deduced by Ph. Rosenau formally, who proposed 
three alternative ways to introduce them [21]. Also, this equation was formally derived by Brenier [T3] 
by means of Monge-Kantorovich's mass transport theory and named relativistic heat equation after him. 
As Brenier pointed out in [M], see also [34], the relativistic heat equation (|f .2|) is one among the various 
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flux limited difl^usion equations used in the theory of radiation hydrodynamics |21j . A general class of 
flux limited diffusion equations and the properties of the relativistic heat equation have been studied in a 
series of papers [SJ IHl El ISl [HI HH] , where the well-posedness of the Cauchy, the Neumann and the Dirichlet 
problem for the relativistic heat equation is proved. 

The above discussion on linear diffusion versus flux-limited diffusion leads to introduce the following 
change in the classical flux 

(1.3) J'=-z/Vu, iy>0, 
associated with the heat equation (or the Fokker-Plank equation) 

Uf = j^Au, 

by a flux that saturates as the gradient becomes unbounded. To do that, it was proposed to link u to the 
flux T through the velocity v defined by the relation T = uv. Along with (|1.3|) this gives 

(1.4) v= ^" 



u 

I Vm I 



According to (|1.4p . if |^| f oo, so will do v. However, the inertia eflFects impose a macroscopic upper 
bound on the allowed free speed, namely, the acoustic speed or light speed c. With this aim, Rosenau 
proposed to modify (|1.4p by taking 

(1.5) V = . 

u ,/^_ Mi 



The postulate (jl.Sp forces v to stay in the subsonic regime (in the case c is the acoustic speed). The sonic 
limit is approached only if | — | ^ oo. Solving (|1.5|) for v, we obtain 

T = uv = 



-1 I / z^lVu 



As we have mentioned before, the motivation for studying the system (|l.ip comes from the transport 
of morphogenes in biological systems. This is a classical problem since the pioneering work of Turing 
[35], Meinhard, Wolpert [3S] or Lander [23. Lander focussed the question as a main problem in the 
understanding of the transport of proteins via signalling pathways: Do morphogen gradients arise by 
diffusion? The relevance of our study is founded on the analysis of the Hedgehog (Hh) signaling which has 
been found to play multiple roles in development, homeostasis and disease (reviewed in [22]). In vertebrates 
the Hh family comprises three proteins (Sonic, Desert and Indian), which act as secreted, intercellular 
factors that affect cell fate, differentiation, survival, and proliferation in the developing embryo and in 
many organs at one time or another. Sonic Hedgehog (Shh) signaling has also an important role in tumor 
formation: the deregulation of the Shh pathway leads to the development of various tumors, including 
those in skin, prostate and brain [25l [26l [31] . The idea is to analyze the morphogenetic patterning of the 
vertebrate embryonic neural tube along the dorsoventral (D-V) axis. The transport of the morphogen Shh 
along the D-V axis in the neural tube represents a natural privileged direction for the description of Shh 
propagation. Actually, the system is symmetric with respect to this axis and this justifies the reduction to 
one dimension. The discussion concerning whether the gradient formation of morphogenes is produced or 
not by diffusion is a central and classic topic in developmental biology. This gives a continuous feedback 
between mathematical modeling and biological experiments, see [20l [28l |35l [19] . Recent results in biology 
provide some findings that really call into question the hypothesis of diffussion which has been so often 
used to model these phenomena: 1) Concerning the cellular differentiation, the role of the quantity of 
morphogen received is as least as relevant as the time of exposure. With linear diffusion models every 
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point (cell) of the neural tube receives instantaneously the information of the niorphogen, [17l [28l [23] . 2) 
Morphogenes are transported in aggregates of several molecules that also include other morphogenes or 
molecules. Then, the typical size of the cluster aggregates is big (of order 1/10) in comparison with the 
extracellular matrix where they are moving [33] . Also their concentration is quite dilute |16[ I33[ I17| . 
Therefore, Brownian motion docs not seem to be the more appropriate choice. 3) In some cases, such as 
with the Hh niorphogen, it has been proved that in absence of another cell-surface protein, called Ihog, 
there is neither propagation nor gradient function of Hh J16j . 4) There do exist privileged ways/paths of 
propagation in the extracellular matrix, a fact that makes the system resemble a traffic map, more than a 
linear diffusion system [T31 [TB] . 

In this setting, the present paper tries to give some insight on this biological problem where the model 
here studied is a first step towards a complete model consisting in 

— -7 — = a{u{t, x), u{t, x)x)j: - f{t - r, u{t, x)) u{t, x) + g{t, u{t, x)), 

where / stands for the concentration of transmembrane receptor in the cells, g represents the concentration 
of the complex binding the niorphogen to the receptor, and where the dependence on u is given through a 
coupling with a system of seven ODE's modeling the rates of change of the concentrations of the proteins 
participating in the signaling pathway coming from the biochemical cascade inside the cells, sec [27j . In 
that work it was also proved that numerical evidence fully agrees with the experiments from a quantitative 
as well as qualitative (propagation of fronts instead of linear diffusion behaviour) point of view, see |161 I30j . 

In addition to the biological or physical motivations, the mathematical analysis of this equation poses 
several difficulties, making even more interesting its study, such as the existence and evolution of fronts 
as well as the study of its finite speed of propagation, the related lack of regularity and the set-up of 
an appropriate functional framework to give a meaning to the differential operator and the boundary 
conditions. In fact, this fiux-limited equation provides a behaviour more related to hyperbolic systems 
than to usual diffusive (Fokker-Plack) systems. To deal with these mathematical problems we need to 
combine and extend the applicability of different techniques stemming from parabolic and hyperbolic 
contexts such as Crandall-Liggett's theorem, Minty-Browder's technique, the concept of entropy solution, 
and the method of doubling variables due to S. Kruzhkov. 

This paper deals with a preliminary study of p.l|) consisting in the analysis of the following system 
(1.6) 



3u 

— = {a{u,u,))^, in ]0,T[x]0,i[, 



-a(u(i,0),M^(i,0)) =/3 >0 and w(t,L)=0, on te]0,r[, 

u{0,x) = uq{x), in a;e]0, L[, 

where the boundary conditions must be considered in a weak sense. Our main result is 

Theorem 1.1. For any initial datum < ito G L°°{]0,L[), there exists a unique bounded entropy solution 
u of il.6\) in Qt =]0, T[x]0, i[ /or every T > such that u{0) = Uq. 

The paper is structured as follows. In the next section we introduce all the tools needed to develop the 
theory: a suitable integration by parts formula, lower semi-continuity results and a functional calculus, in 
order to be able to give a sense to the differential operator. In Section 3 we discuss the associated elliptic 
problem: we define what a solution is, and then we prove existence and uniqueness of such a solution. 
Next, this material is used to define an accretive operator and construct a nonlinear semigroup, which 
accounts for solving (|1.6p in a mild sense; all this is the content of Section 4. In Section 5 we prove 
that the mild solution previously constructed can be characterized in more operative terms, as a so-called 
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entropy solution -a concept which is also introduced in this section-, and we prove a comparison criterium 
which in particular entails uniqueness of entropy solutions, thus proving Theorem II .11 

2. Preliminaries 

2.1. BV functions and integration by parts. For bounded variation function of one variable we follow 
[T]. Let / C M an interval, we say that a function u € L^{I) is of bounded variation if its distributional 
derivative Du is a Radon measure on I with bounded total variation \Du\{I) < +oo. We denote by BV{I) 
the space of all function of bounded variation in /. It is well know (see [1]) that given u £ BV{I) there 
exists u in the equivalence class of u, called a good representative of u with the following properties. If J^ 
is the set of atoms of Du, i.e., x £ Ju if and only if Du{{x}) ^ 0, then u is continuous in / \ J^ and has a 
jump discontinuity at any point of J„: 

u{x-) := limu(y) ~ Du{]0,x[), u{x+) := limu(y) = Du{]0,x]) \f x E J„, 

ytx ylx 

where by simplicity we are assuming that / =]0, L[. Consequently, 

u{x+) — u{x-) = Du{{x}) y X S Ju- 

Moreover, u is differentiable at C^ a.e. point of /, and the derivative u' is the density of Du with respect to 
£^, being C^ the d-dimensional Lebesgue's measure. For u G BV{I), the measure Du decomposes into its 
absolutely continuous and singular parts Du = D°^'^u + D-''u. Then Z3°'^u = u' C^ . Obviously, if u G BV{I) 
then u S M^^'^(/) if and only if D^u = 0, and in this case we have Du = u' C^. From now on, when we 
deal with pointwise valued i?F-functions we always shall use the good representative. Hence, in the case 
u £ W^'^{I), we shall assume that u £ C{I). 

Given z S iy^'^(/) and u E BV{I), by zDu we mean the Radon measure in / defined as 

{ip,zDu) -.^ ifizDu V^e Cc(]0,i[). 
Jo 

We need the following integration by parts formula, which can be proved using a suitable regularization 

of u S BV{I) as in the proof of Theorem C.9. in [5]. 

Lemma 2.1. If z e W^-^{I) and u e BV{I), then 

zDu+ / u{x)z'{x)dx = z{L)u{L^)-z{0)u{0+). 
Jo 

2.2. Properties of the Lagrangian. Hereafter C denotes a generic constant, its value may change from 
line to line. We define 

(2.7) .i..,0:^^M=. 

We assume a(z,0) = for all z G K. Then a(z,^) — d^F{z,£^), being the Lagrangian 



By the convexity of F, 

(2.8) aL{z,^){r]-^)<F{z,r])-F{z,^) for aU z,^, ?? e R 
Note that we have 

(2.9) c\z\\^\ z^ <a{z,0^<cM\^\ foranz,^eM, \z\ < M. 

V 
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Moreover, using ()2.8p it is easy to see that 

(2.10) (a(z,C)-a(z,^")).(C-e)>0 

for any (z, ^), (z, C) € M x M, \z\ < M. 

We introduce the following notation to ease the way in which our functional calculus is written: for any 
function q let Jq{r) denote its primite, i.e., Jq{r) = L q{s)ds. 

Assume that /:]RxM^-[0,cx)[isa continuous function convex in its last variable such that 

(2.11) < f{z, < C(l + ICI) V(2, e R X M, \z\ < M. 

for some constant C > which may depend on M. Given /(z, ^), we define its recession function as 

r 



r(z,0- lim tf I z 

t^Q+ \ t 

We assume that /"(z,^) = Lp{z)ip'^{£^), with ip Lipschitz continuous, f/;" homogeneous of degree 1. Then, 
working as in |S], if for a fixed function G Cc{]0,L[) wc define the operator 7?,^/ : BV{]0,L[) ^- M by 

(2.12) n^f{u):= j^ c^{x)f{u{x),u'{x))dx + J^ ^(^:r.)^jj" (^\ \D' ^u)], 

wc have that TZ^f is lower semi-continuous respect to the L^-convergence. 

For instance, wc discuss here for future usage one of the most recurrent cases: define 6[z) — c\z\, and 
note that F°{z,^) = 0{z)'iP'^{^), with V°(0 = l?l- Therefore, 

n^Fiu):= [ ^{x)F{u{x),u'{x))dx + ^ j ^{x)\D'{u^)\ 
Jo ^ Jo 

is lower semi-continuous in i?y(]0,L[) respect to the L^-convergence. 

We shall consider the function /i : M x M ^ M defined by h{z, ^) := a(z, ^) • ^. Note that 

(2.13) h{z,^)>0 VC,zeM. 
We will make use of the following property: 

(2.14) h°{z,^)=F°{z,0 Ve,zGK. 

As for the Dirichlet problem (see ^U\), in general, the data in L is not taken pointwise; we need to 
introduce functional that take into account the boundary. The following result is a particular case of 
Theorem 2.4 in [TU] 

Theorem 2.2. Let f be verifying \2.11\) and /"^(z,^) = 99(z)|^|, with ip Lipschitz continuous, let (/) S 
C([0,i])+ be given. Then, the functional J-'^f : BV{]0, L[) — > M defined by 

J^°j{u) -.^n^fiu) + c^{L)\Mu){L^)\ 

is lower semi- continuous with respect to the L^ — convergence. 
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2.3. Spaces of truncated functions and associated calculus. We need to consider the following 
truncature functions. For a < b, let Tafi{r) := max(niin(6, r), a). As usual, we denote Tk = T-k.k- We also 
consider the truncature functions T^ j,(r) := Tafi{r) — I {I E R). We denote 

Tr := {Ta,b : < a < 6}, T+ := {T^^i, : < a < &, ? € R, T^^^ > 0}. 
Given any truncature function Tfc, we define 

Tkir)+ := max{Tfe(r),0} and Tkir)- := min{rfe(r),0} = -Tfc(-r)+, r e R. 

Consider the function space 

TBV+{I):^{ueL\l)+ : T{u) e BV{I), y T e %} ; 

we want to give a sense to the Radon- Nikodym derivative u' of a function u e TBV'^{I). Using chain's 
rule for BV-functions (see, for instance, [T]), and with a similar proof to the one given in Lemma 2.1 of 
PT| . we obtain the following result. 

Lemma 2.3. For every u G TBV^{I) there exists a unique measurable Junction w : / — > R such that 

(2.15) {Ta,b{u))' = vX^a<u<b\ r'-a.e., V T^.b e T.. 

Thanks to this result we define u' for a function u G TBV^^) as the unique function v which satisfies 
(|2.15p . This notation will be used throughout in the sequel. The notation dx will also be used in the 
case of functions of several variables (say t and x), for the same purposes, whenever there is some risk of 
confusion. 

We denote by V the set of Lipschitz continuous function p : [0, +oo[— > R satisfying p'{s) = for s large 
enough, and write V^ := {p E V : p> 0}. We recall the following result ([5], Lemma 2). 

Lemma 2.4. If u E TBV^{I), then p{u) E BV{I) for every p G V such that there exists a > with 
p{r) ~ for all < r < a. Moreover, with the above notation [p{u)]' = p'(u)u' C^-a.e. 

For u e TBV+{]0,L[) we will define 

m(0+) := lim Ti„(w)(0+) and w(L_) := lim Ti „(m)(L_). 

It is easy to see that the above limits exist. 

Let S" e 7^+ andT = T^,,. Given w G TBV+{]0,L[), Lemma[2l]assures that S{u)T{u), Jt's{u), Jts'{u) G 
BV{]0,L[). Moreover, D{S{u)T{u)) = DJt's{u) + DJts'{u) and hence, if z G W^^^{]0,L[), 

zD{T{u)S{u)) = zDJx'siu) + zDJrs'iu). 

FoiueTBV+{]0,L[), (f>e Cc(]0,i[), T ^Ta^-le r+ and / as in the previous subsection -see dSIIII)-, 
we define the functional 



'R{cj>f,T){u):^n^f{Ta,b{u))+ / ct>{x){f{u{x),Q) ~ f{a,Q))dx 

[ii<a] 



^(a;)(/(w(x),0)-/(6,0))da;. 

[u>b] 

We have that TZ{4if,T){-) is lower-semi-continuous in TBV^{]0,L[) with respect to the i^-convergence. 
Given S,T £ T^ and u G TBV^{\{), L[), we define the following Radon measures in ]0, L[, 

{F{u,DT{u)),cf>) :=n{^F,T){u), {Fs{u,DT{u)),cf>) :^n{4>SF,T){u), 
{h{u, DT(u)),4>) := n{4>K T)iu), {hs{u, DT{u)), 4>) := 7^(05/l, T){u), 
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for (/) e Cc(]0,i[). Using ([2l^ and (f234l) . we compute 

F{u,DT{u)y = ^ \D'{T{u)f\ = h{u,DT{u)y, Fs{u,DT{u)y = \D'Jse{T{u))\ = /i5(^, ^n«))^ 

;i(u, DTiujf = h{u, (r(u))'), /^s(it, DT{u)Y'' = s-CiO/iCu, (r(u))'). 

3. The Elliptic Problem 

Given v e i^(]0, L[), we are interested in the following problem: 

[ -(a(w,u'))' =v in ]0,L[ 

(3.16) <^ 

[ -a(M, u')U=o = /3 > and u(L) = 0, 

where a is given by (|2.7p . We introduce the following concept of solution for problem p.l6p . 

Definition 3.1. Given w € L^(]0,L[), we say that u > is an entropy solution oi <l(i.lQ\ if it e ri?F+(]0, L[ 
and a(M, u') G C([0, L]) both satisfy 

i> = -7:>a(u,M') in X''(]0,L[), 

— a(w, i/)(0) = /3, and a(M, u'){L) = —cu{L^). 

(3.17) /i(w, i:'r(M)) < a(w, u')DT{u) as measures VT e r+ 



+ 



(3.18) hs{u,DT{u))<a{u,u')DJT's{u) as measures V5 G 7^+, T e T 



Note that p.l7p can be rewritten as h{u,DT{u)Y < [a.{u,u')DT{u)Y, and thus it is equivalent to 

^\D%{T{u)y)\ < a.{u,u')D'T{u) as measures VT e r+. 

Also we have that p.lSp can be rewritten as hs{u,DT{u)Y < [a{u,u')DJT'siu)Y , and is equivalent to 

\D%Js0iT{u))\<SL{u,u')D''JT'siu) as measures ySeV+,TeT+. 
Observe that since — a(u, u'){0) ~ (3, we have 

(3.19) m(0+) > - > 0. 

c 

We introduce now the main result of this section. 

Theorem 3.2. For any < f G L°°{]0,L[) there exists a unique entropy solution u £ TBV'^{]Q,L\) of 
the problem 

{u~ {a{u,u'j) =/ in]0,L[ 

-a(u,u')U=o = /?>0 u{L) = 0, 
which satisfies ||w||oo < -^^(Z?) c, z^, ||/||oo)- 
Moreover, let u,u be two entropy solutions of 113. 20]) associated to f,f G i^(]0, £[)"'', respectively. Then, 

L i-L 

{u - u)+ dx< if- /)+ dx. 

"'0 
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Proof. Existence of entropy solutions. Wc divide the existence proof in different steps. 

Step 1. Approximation and basic estimates. Let < / G L°°{]Q,L\). For every n £ N, consider a„(z,^) 

:— a(z,^) H — ^. As a consequence of the resuhs about pseudo-monotone operators in J15j we know that 

n 
Vn e N there exists a unique m„ e Ty-'^'^(]0, L[) such that m„(L) = and 

pL pL -I pL pL 

(3.21) / Unvdx+ / a.{un,u'j^)v' dx ~\ — / wjjw'dx — /3i)(0) = / fvdx 

Jo Jo "■ Jo Jo 

for ah V e M^i'2(]0,i[), w(L) == 0. 

The fohowing resuh can be easily obtained by multiphcation by u~ and integration over [0, L]. 
Lemma 3.3. The functions Un are non-negative Vn G N. 

Now we give a bound for the sequence u„ at zero. 
Lemma 3.4. The sequence {wn(0)} is bounded. More precisely, 



< u„(0) < < 






c V c 



if c> \fv 
if c < \fv. 



Proof. Taking v — u„ in p.2ip . we get 



(3.22) 



1 



+ a(u„,u^Ju^ + -((u„)') (ix = ^u„(0)+/ /it„ dx 



Then, dropping non-negative terms and using Young's inequahty, we get 
(3.23) / u\dx<\ f^dx + 2/3un{0). 



Now we can write u„|u^| = ^Ku^)'], and taking into account (|2.9p we have u^a(M„,M^) > cm„|u^| — ^i 
Then, from p.22p . we obtain 



(3.24) 



l\iuiy\ + ^-^^^^)dx< 

v2 "■ / JO 



- 1 ) u^j + /u„ ) dx + /3Un{0) 



Assuming now that 1 > 0, we apply Young's inequality in the right hand side of p.24[) . which now 

reads 

2 



u: 



^ dx 



f'dx + ^UniO). 



As c > ^/l' we have 77 ^ 5 > 0, which allows us to bring in (|3.23p . thus obtaining 



(3.25) 

Then, we have 



^ \{ulY\dx+- {{ujfdx<- /2d.T+^^«„(0). 







l\uim = l\ui{L)-uim^l 



ujdx 







< 



f^dx 



2l3c^ 



w«(0). 
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from where we get that ^^(O) - ^u„(0) - ^\\I\\l < 0. Hence, for all n e N, 

In case that c^/iy — 1 < 0, from p.24p we obtain 

^ / \{uly\dx + - f {{u^)'fdx< [ fundx + ^UniO). 
2 Jo nJo Jq 

Then, using Young's inequality and having in mind (j3.23p . we get 

(3.26) ^ [ \{ulY\dx + - [ {{uny)^dx< f fdx + 2(3un{0). 

Thus, we have that w^(0) — ^/3u„(0) — f /q /^ < 0, from where it follows that for aU n G N, 



»-"<»> 4 (?W(iy->i)^?-/>^- 



D 



By (inSSl), ([5:^ and Lemma [3H we get 



(3.27) I I \{ul)'\dx+- I {{un)'fdx<C VneN. 

2 Jo '^ Jo 

Lemma 3.5. The sequence {u„ : n G N} is uniformly hounded in L°°{Q,L). 

Proof. By Lemma [3.4[ we know that AI = max{||/||oo,max{u„(0) : 7i E N}} is finite. Then, taking 
V = {un — AI) as test function in p.2ip . it is easy to see that ||u„||oo < M and Lemma [23] holds. □ 

Lemma 3.6. The sequence {u„} is uniformly bounded in TBV'^QO, L[). Furthermore, there exists a 
function < m G TBV^QO, L[) D L°°(\0,L[) such that (up to subsequence) Un ^- u a.e. and strongly in 
LH]0,L[). 

Proof. By Lemma 13.51 extracting a subsequence if necessary, we may assume that u„ converges weakly 
in L^(]0,L[) to some non-negative function u as n — > +00. Moreover, by Lemma [5751 again, we have that 
< u G L°°(]0, L[). On the other hand, if < a < 6, by the coarea formula and p.27p . we have 

\{Ta,b{u^))'\dx^ |i?X[„„<,]|(]0,L[)dt= / |i?X[„2<t2]|(]0,L[)di 

J a J a 

|TO[„.<,,]|(]0,iD:^<;^ / \{ulY\dx<-. 
' "■- " 2^/s 2a Jo a 

Consequently, we may assume that m„ converges almost everywhere to u. Then, by the Vitali Convergence 
Theorem, we get that u„ — > m in L^{]0, L[), and using the above estimate on the gradients we obtain that 

ueTBV+{]0,L[). a 

Since |a(u„,M^)| < c|w„|, by Lemma 13.51 we may assume that 

(3.28) a{un,u'n)^z as n ^ cx), weakly* in i°°(]0, L[). 
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By assumption we have that a(w„, mJJ = c|u„|b(u„, wJJ with |b(M„, u'^)\ < 1 (independent of n), ||it„||oo < 
M, and Un — > u a.e. as n — > cxd, so we may assume that b(u„, u^) ^ Zh as n — ^ oo, weakly* in L°°(]0, L[), 
and 

(3.29) z = cuzb, with ||z6||oo < 1- 
On the other hand, by ((X?f)) . 

(3.30) -<^0 in L^{]0,L[). 

n 

Given (p E 2?(]0,i[), taking w = ^ in p.2ip we obtain 

pL -I pL pL 

Un(f>dx+ / a(w„, w^j)0' dx- H — / u'^(j)' dx ~ / f(j)dx 
^0 "Jo Jo 

Letting n — > +oo, having in mind p.28p and p.30p . we obtain 

/ (/ — u)(j) dx — z ■ (j)' dx, 
Jo Jo 

that is, 

(3.31) f-u=-Dz, in I?'(]0,L[) 

and 

(a„(u„,0)' ^Dz weakly in L^(]{),L[). 

Note that by (P3T|) . we have z e Tyi'i(]0, L[) and Dz = z'. 

Working as in the proof of Lemma 4.2 of [5], we can prove the identification 

(3.32) z{x)=Si{u{x),u'{x)) a.e. xe]0,L[ 
From ([332]) and (f33T1) it follows that 

f -u = -D&{u,u'), inX>'(]0,L[) 
Lemma 3.7. T/ie /Z-ux — a(w, m') verifies the Neumann condition at x = 0. 

Proof. Let w e W^'^{]0,L[) such that w{L) = and consider Wk € W'i'2(]0,L[) with WkiL) = for ah 
fc € N, lUfc — > w pointwise and wj. — )■ w' in L^(]0, L[). Taking in (|3.2ip Wk as test function and letting 

n — > +00, we get 

pL pL pL 

/ uwkdx+ / zw'f^dx — (iwk{Q) ^ / fwkdx. 
Jo Jo Jo 

Then, letting k — > +00 we arrive to 

(3.33) / uwdx+ / zw' dx — I3w{0) = / fwdx. 

Jo Jo Jo 

Fixed w e By(]0,L[) such that w{L_) = 0, let w,„ e M^i'i(]0,L[) with ix;„,(L) = 0, w;,„(0) = w(0+), 
and such that Wm — > w in L^(]0, i[). Taking in p.33p Wm as test functions and integrating by parts we get 

(^L pL pL 

{f — u)w,ndx~ / zw'jj^dx — f3w{0+) = — / z'wmdx — w{0+){z{0) + (3), 
Jo Jo 

and letting m — > +00, we obtain — z(0) = /3. □ 
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Lemma 3.8. Let S eV+,T eT+ and (j) e C\[0,L]), > 0, with (f>{0) = 0. Then 

'^ c 

ct^F{u,DT{u)) + ct^{L)-\{T{u)f{L^)\ 



(3.34) 



and 



(3.35) 



In particular, 



<l <PzDT{u)+ (j)F{u,0)dx-(l){L)T{u){L^) + <j>{L)\Je{T{0))\ 



L 

(j>Fs{u,DT{u)) + ^{L)\Jgs{T{u){L_))\ 



< / (I>zDJt's{u)+ / (f>S{u)F{u,0)dx 
Jo Jo 

-4>{L)z{L)Jt's{u{L-)) + 0(L) \Jgs{T{{)))\ . 



(3.36) F{u, DT{u)) < zDT{u) + F{u, 0)C^ as measures in ]0, L[. 

(3.37) Fs{u,DT{u)) < zD{Jt's{u))) + S{u)F{u,0) £^ as measures in ]0,i[. 

Proof. We will only prove dSSS]), the proof of (|XM)) being similar. Let < e CH[0, L]) with 0(0) = 0. 
Since -F^gp is l.s.c. with respect to the L^-convergence, letting n — >■ oo we obtain 

cf>Fs{u,DT{u)) + ^{L)\Jes{T{u){L^j)\ 





<liminf / cj)Siu„)Fiu„,T{uJ) dx + ^{L) \Jgs{T{0))\ 
<limsup / cj)Siu„)F{u„,T{u„y) dx + <j){L)\Jgs{T{0))\ 

Tl— !-00 JO 



By the convexity ()2.8|) of F and using that a{umT{un)')T{un)' = a(u„,u^)T(w„)', we have 



L 

(pS{Un)F(Un, TiUn)') dx 


L 



< / 4>S(Un)ai{Un,T{Un)')T{Un)' dx+ j (pS {Un)F{Un,0)dx 

Jo Jo 

<f>aL{Un,u'^){JT's{Un)y dx + / (/)S{Un)F{Un,0)dx. 



Now we take v = JT'si'an)4' ^s test function in (j3.2ip ) and we obtain 



1 '•^ 



0a(u„,M^)(JT'S'(Mn))'rfa;+ - / (J)u'^{Jt' siun))' dx 



"^ Jo 

L rL -| pL 

{f -Un)JT's{un)(pdx- JT's(wn)a(w„,w'J(/)'da; / JT's{un)u'^4>' dx. 

Jo "IT- Jo 



D 
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Letting n — > cx) we get 

limsup / (j)a.{un,u'n){JT's{un)y dx < (j){f - u)JT's{u)dx - JT's{u)z(t)' dx 

n Jo Jo Jo 

= / (bzD{JT'siu)) - (l)iL)z{L)JT'sHL_)). 
Jo 

Finally, 

/ <l>Fs{u,DT{u)) + <j>{L)\Jes{T{u)){L^)\< f ^zDJt's{u) 
Jo Jo 

+ (/.(L)|Jes(T(0))|-(/)(L)z(L)JT's(w(i-))+ / (l)S{u)F{u,Q)dx 

Jo 

and (p:^ holds. 

Lemma 3.9. The inequalities {3. 17^ and hS.l^) hold. 

Proof. Using (|3.36p and the fact that h{u,DT{u)) is a measure concentrated in ]0,i[, it follows that 

h{u,DT{u)y = F{u,DT{u)y < (z DT{u)y. 

Hence, 

zDT{u)) = zT{u)'C^ + (zDT{u)Y > zT{u)' C^ + h{u, DT{u)y = h{u, DT{u)), 
and (|3T7l) holds. 

Using p.37p we have 

zD{Jt's{u))) = {zD{JT's{u))y^ + {zD{JT's{u))y > z{Jt's{u))' + {Fs{u, DT{u))y 

= z{JT's{u)yC^ + {hs{u, DT{u))y = hs{u, DT{il)), 

and wc obtain (|3.18p . □ 

Lemma 3.10. The Dirichlet condition a(u,u')(L) ~ —cu{LJ) holds. 

Proof. Firstly, observe that by (|3.29p we have 

\z{L)\ < cu{L^). 

Then, it is enough to prove the lemma in the case u{LJ) > 0. In that case, again by p.29p and having in 
mind that z is continuous in [0, L], we have 

(3.38) z(L) =cw(L_)C, with |C| < 1- 

Given T G T+, for tti > 1 wc consider S := y"-i <= 'p+ ^ Taking singular parts in p.35p wc have 

(3.39) \JeT^-i{T{uy{L-)\ < -z{L)Jt^-it'{u{L^)) + \JeT'^.-i{T{0))\ . 

Consider now T = Td,d' with < d < u{L-) < \\u - ||oo < d'. Using (|338l) . the inequality (|339l) 
particularizes to 

£^"+1 + —^— (u^+Vi^) - d^+i) < -d"'+^ - -Cu{L) [vJ^iL^ - d") 
2 TO + 1 2 m 

and letting d ^- 0+ wc have 



"= u^+\L^) < --u{{L^))iu^{{L^). 



m+1 
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Then, since u{L-) > 0, we get ;j^ < — ^ for all 1 < m. Therefore, since |^| < 1, we have ^ = —1. 
Consequently, by p.38p we finish the proof. □ 

Proof of uniqueness. Let u,u be entropy solutions of (j3.20p associated with /, / e L^(]0, L[)^; respectively. 
Let pn be a classical moUifier in ]0, i[, ij] G 2?(]0, L[) and & > a > 2e > 0. Let us write 

£.n{x,y)= Pn{x-y)i^{-^^\ , and r = T^,,. 

We need to consider truncature functions of the form Se,i{r) :~ T^{r — /)+ = Ti^i^^{r) — I e T^ and 
S[{r) ■= T^{r ^l)- +6 = Ti^e.iir) + e - I e T+ , where I > 0. Observe that Si{r) = -T,{l - r)+ + e. 

If we denote z(y) = a.{u{y) , dyu{y)) and z{x) = a.{u{x) , dxu{x)) , we have 

w-z' = / and u-z =J in V'{]0,L[). 

Then, multiplying the equation for u hyT{u{y))Sf^u(^x){u{y))£,nix,y), that ioi u hy T {u{x)) Se {u{x))£,n{x , y), 
integrating in both variables 

/ [u{y)T{u{y)) - u{x)T{u{x))]T^{u{y) ~ w(a;))+^„(a;, y) dxdy 
Jo 

el / {u{x) - f{x))T{u{x)) + ^„(x, y) dxdy 
Jo 

L .L 

/ Ux,y) (zDy[T{u)S,^j,^,jiu)]dx + zDx[T{u)Sfy'>{u)]dy 
Jo ^ 

L pL 

/ T{u{y))S^^u(x)iu{y))z{y) ■ dy^n{x , y) dxdy 
Jo 

+ / T{u{x))Sfy\u{x))zix)-d,Ux,y)dxdy 

Jo Jo 

L nL 

[f{y)Tiuiy)) - J{x)T{u{x))]T4u{y) - u{x))+^n dxdy. 

Let / denote all the terms at the left hand side of the above identity, but the first one. From now on, since 
u, z are always functions of y, and m, z are always functions of x, to make our expressions shorter, we shall 
omit the arguments except in some cases where we find useful to remind them. 

With slight modifications of the method used in the proof of uniqueness in [5] we can obtain the following 
result. 

Lemma 3.11. The following inequality is satisfied 
e 



(3.40) 



I>o{e)- i I / ir,zD,T{u)\dy+-j J T,{u ~u)+{T{u)z - T{u)z) ■ {d,in + dy^,,) dx dy 
where o(e) denotes an expression such that o{e) — > as e — > 0. 
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By the above lemma, dividing p.40|) by e and letting e — ?► we obtain 

f-L /-L 

£.n{x,y){u{y)T{u{y)) -u{x)T{u{x)))sign+{u{y) -u{x))dxdy 

"'0 

/ Pn{x~y)s\gTi^{u{y)-u{x)){T{u{y))z,{y)-T{u{x))z{x))-^/ {-^^\ dxdy 

pL pL pL I pL \ 

</ / Ux,y){f{y)T{u{y))~J[x)T{u{x)))s\gn+{u{y)-u{x))dxdy+ / Ux,y)^D,T{u)\ dy, 

Jo Jo Jo \Jo J 

where 

. +, . ( 1 ifr>0 

Letting n — s^ oo, we find 

/ ip{x){u{x)T{u{x)) — u{x)T{u{x)))signQ {u{x) —u{x))dx 
Jo 

signQ (li(a-) — u{x)){T{u{x))z{x) — T{u{x))z{x)) ■ -0' {x) dx 

< I ij{x)[f{x)T{u{x))-J{x)T{u{x))]sign^{u{x)~u{x))dx+ I ip{x)2.DT{u). 
Jo Jo 

Taking now a sequence ipm t ^]o.l[, i^m G T^Q^^Ll) in the above formula, we have 

{u{x)T {u{x)) — u(x)T(M(x)))signJ(u(a:) — ^(a:^)) dx 
'a 

+ fim / aign^ {u{x) ~u{x))(T {u{x))z{x) — T{u{x))z{x)) ■ ijj'^^x) dx 

pL pL 

< {f{x)T{u{x))~7{x)T{u{x)))s\gn+{u{x)^u{x))dx+ zDT{u). 
Jo Jo 

Now we deal with the second term in the above expression. 

lim / sign(t(u(a;) — u{x)){T{u{x))z{x) — T{u{x))z{x)) ■ ^jj'^^ (x) dx 

= — lim / ■ipm{x){zD[signQ{u~u)T{u)]~z{x)D[s\gnQ{u~u)Tiu)\\ 

+ lim / V'm(a;){sign^(M(x) — u{x))T {u{x))z' (x) — sign(t(?i(a;) — u{x))T(u(x))z' (x) } dx , 
which leads to 

signQ {u{x) — u[x))T {u{x))z' [x) dx — I zD[sigi\Q {u — u)T {u)\ 
Jo 

L i-L 

signj(w(x) — u{x))T {jl{x))z' {x) dx + zD[signQ (u — u)T{u)] 
Jo 

= [z{0)T{u{0+)) -ziO)Tiu{0+))]sign+{u{0+) -u{0+)) 
[z{L)T{u{L^)) - z[L)T{u{L^))] sign+(w(i_) - u{L^)). 
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Therefore, 

L 

{u[x)T {il{x)) — u{x)T {u{x)))signQ [u[x) — u{x)) dx 



+ [z(0)r(u(0+)) - z(0)r(ll(0+))] sign+(i.(0+) - 11(0+)) 
- [z(L)r(u(L_)) - z(L)T(^l(L_))] sign+(^(L_) - u{L^)) 

pL pL 

< [f{x)T{u{x))-J{x)T{u{x))]sign+{u{x)-u{x))dx+ zDT{u). 
Jo Jo 

Dividing by 6 > 0, and letting a — > 0^, and letting 6 — > 0^ in this order, we obtain 
[uXiuyo] - u^[u>o])sign([(u - u) dx 

+ [z(0)sign+(7.(0+)) -z(0)sign+(il(0+))] sign+(w(0+) - 11(0+)) 
- [z(L)sign+(u(L_)) -z(L)sign+(ll(L_))] sign+ {u{L^) -u{L_)) 

< / (/^[«>o] - f^[u>o])signo {u - u) dx + lim - lini / zDT{u) . 
Jo b^o b ya^ojg J 

Now, since z(0) = z(0) = -(5^0, and u(0+) > | > and u(0+) > | > by (plQ)) . we have 
that the second term in the above expression vanishes. On the other hand, since z(L) = —cu(L-) and 
z(L) = —cu(L^), the third term in the above expression is non-negative. Consequently, 

("^[u>o] - MX[jj>o])sign([(w - u) dx 

JO 

(3.41) < J (/X[„>o] -7%>o])sign+(u-IZ)da;+ hm i Mimy zDT{u)] . 

Next we claim that 

(3.42) / = a.e. on [u = 0] and 7 = a.e on [m = 0]. 

Let < e ^CO,i[) be and a > 0, e > 0. Multiplying / - u = -z' in V'{]0,L[) by T;^,a+eiu)<t> and 
integrating by parts and having in mind p.l7p and (j2.13p . we have 

(/ - u)T^\,+,{u)c^dx = f c^zDT^\,^,{u) + / z • 0'T,%+,(u) dx > f z • 0'r,%+,(u) dx. 
'0 Jo Jo Jo 

Dividing by e and letting e — >■ 0+, we get 

(/-'u)X[„>a]0dx > / z • ^'X[„>a] da;. 
JO Jo 

Hence 

fL pL pL 

if -u)^[u<a](l)dx = {f-u)(j)dx- [f -u)X^^ya]{x)(l)dx 
Jo Jo 

pL pL pL 

< / if -u)(t>dx- z • (/)'X[„>a] da- = / z • </)'X[„<a] dx. 
Jo Jo Jo 
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Then, letting a — > 0+, since z = in [u = 0], we have 

PL, pL 

/ .f^[u=Q]<pdx = (/ - u)X[^^o](l)dx < 0, 

Jo JQ 

for all < G T-^Q^j^l), from where it follows that /^[u=o] = a-C- in ]0,L[. Similarly, fX^j-Q^ — a.e. 
in]0,i[and (lOS]) holds. 

On the other hand, by (|3.42p we have, 

lim - I lim / zDT(u) ] 

b^O b \a^0j^ ' j 

= - lim i lim I z(0)T(w(0+)) - z{L)T{u{L^)) + [ T{u)z' dx \ 

b^O b a^O y Jo J 

= - lim i I z(0)ro,fc(u(0+)) - z(L)ro.faNi-)) + f n,b{u)z' dx J 

= -z(0)sign|]"(?l(0+)) + z(L)sign(]"(u(L_)) - / X[:ej>o]z' dx 

Jo 



--z(0)sign+(u(0+)) + z(L)sign+(ll(L_))- / z' dx 

JQ 

= z(0) (1 - sign+(^(0+))) + z(L) (sign+(il(L_)) - l) = 0. 
Then, from (|OT|) . it follows that 

/ iuX[u>o] - uX[^j>o])signJ(u ~u)dx < / (/X[„>o] - /^[«>o])sign|;!"(w - u) dx. 
Jo Jo 

Hence, using p.42p . we obtain 

pL pL pL 

/ {u-u)+dx< / {f-J)sig4{u-u)dx< / if-D+dx. 
Jo Jo Jo 

This concludes the proof of the uniqueness part of Theorem 13.21 D 

4. Semigroup solution 

In this section we shall associate an accretive operator in i^(]0, L[) to the problem (|3.16p . 

Definition 4.1. {u,v) e Bp if and only if < u e TBV^(]0,L\), v e i^(]0,i[) and u is the entropy 
solution of problem p.l6p . 

From Theorem l3.2[ it follows that the operator Bp is T-accretive in L^{]0,L[) and verifies 

(4.43) L°°{]0,L[)+ cR{I + XBi3) for aU A > 0. 

In order to get an i°°-estimate of the resolvent, we need to find the steady state solution, that is, the 
function up which is the entropy solution of the problem 

f -(a{up,u'^))' ^0 in ]0,L[ 

(4.44) I 

[ -a(ii/3,w^)|^=o = /3 > and w/^(L) = 0. 
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Proposition 4.2. There is a non-increasing function up G C^{]0, L[), with up > -, that is an entropy 
solution of the stationary problem ^.44\ l- Moreover, there exists a constant M :— M(c, (3,1/, L) such that 

\\up\\ao<M. 



Proof. 

Integrating (|4.44|) over ]0, L[ we find that a{up,u'p){L) — ~/3. Now, if up has to fulfil the weak Dirichlct 
condition a{up,u'p){L) ~ ~cup{L-) then we must have up{L-) — (3/c. We will follow this prescription 
hereafter. 

If Up is a solution of the problem (14.441) , we have 



{a{up,up)) =. ^^=^ V 



upup 



fiup 



Then, assuming that m^ < 0, we get 



Thus, we get that up satisfies the ordinary differential equation 



Up\u-^- It 



up 



By means of the change of variable 't^^ = "a ~ ( f ) > we arrive to the ODE 




Then, 



dy^l v'{y) dy 



v'iy) 



■dy 



vjy) 
p/c 



, , /3 fv(L)\ (3 fvixY 

v(L) — v[x) arctan H — arctan — — — , . 

c \l3/c J c Xfi/c J 



Hence, we get 
(4.45) 



X — L — —\\ upixY — — -A — arctan 
pV Vc/ c 



^Jupixr-^^^ 



li X — Ug^(y), then we can write (|4.45p as 



".-<»' ^-ivRF^^ 



arctan 



-^r 
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Thus, 



u-,') {y) 



y2_(l 



■2 \/3j \ * ' c2y2 



-1 



and consequently, since {up){L^) = |-, we obtain that 



1 ,. V^'^(^)V/3\^ cv 



(u^)'(i_) = hm — = hm . , . 



^^' (-^') (y) '^^ 



Finally, since up satisfies — (a(u^(a;), Wo(a;)))' = if a; g]0,L[ and satisfies the boundary conditions also, 
we have that up is an entropy solution of the problem ()4.44p . 

D 

The following homogeneity of the operator Bp will be important to get the L°°-cstimate of the resolvent. 

Proposition 4.3. For yU > 0, A > and /3 > 0, we have 

(4.46) {I + \Bp)-\^iu) = ^l[I + \Bs^ \u). 

Moreover, for (3i < ^2, u ^ L°°{]Q,L[)+ and X > such that {I + XB(}.,y^ (u) e BV{]0,L[), we have 

(4.47) {I + XBf,,y\u) < {I + XBf,,y\u) a.e.xe]0,L[. 

Proof. From the definition of the operator it is easy to see that if w G D{Bfi_), then jiu G D{Bp) and 
Bp{iJ,u) = fj,B^{u). Then, we have 

«:=(/ + XBr) (/iw) <;=^ V + XBr{v) ~ fiu <(=^ —V H — XB Av) = u 



<^=^ -v + XBji(-)^u -^^ (l + XBjA (u) = -, 

from where (|4.46p follows. 

Finally, let us see that (|4.47p holds. Let Ui := (/ + XBp.)^ (u), i = 1,2. Then, Ui is an entropy solution 
of the problem 

Mi — A (a(ui,u9) =u in ]0,i[ 

-a(ui, u-)U=o = A > and u{L) = 0. 

Therefore, if p„ are non negative increasing functions that are an approximation of the sigup function, 
having in mind p.lOp . since p„(ui — U2) € BV{]0,L[), we get 

/ {ui - U2)pn{ui - U2) dx = x({a{ui,u[))'-{a{u2,U2))']pn{ui-U2)dx 

= X{a.{ui,u[) - a{u2,U2)) D{pn{ui - U2)) 

Jo 

+A(a(ui,ui)(L_) -a{u2,u'2){Ly)pn{ui - U2){Ly - X{a{ui,u[){0+) - a{u2,U2){0+)) Pn{ui -U2)(0+) 

< X{a{ui,u[){Ly-a{u2,u'2){Ly)pn{ui - U2)(i-) + A(/?i - l32)Pn{ui -U2)(0+) 

< A(a(ui,ui)(L_) -a(u2,U2)(L-))p„(ui -U2)(L_). 
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Then, taking limit as n — >■ +cx) we get 

(mi -U2)^ dx < A(a(ui,Ui)(L_)-a(u2,U2)(-^-))sign(t(wi --U2)(i-) < 0, 

since a(wi,u^)(L) = —cui{L^), i ~ 1,2. Therefore, ui < U2, and wc finish the proof. □ 

Proposition 4.4. For u E L°°(]0,i[)+ and X > 0, we have 

< {I + XB^y^ {u) < iiup, with At = max|^^i^^,l 

Proof. Let up be the entropy solution of the stationary problem (|4.44p given in Proposition 14.21 Then, 
(u/3,0) e Bjs, from where it follows that 

(4.48) {I + XBp)-^ (up) = up. 

On the other hand, since up > §, if /x := max{ "^"^^'"^ ^ 1}^ we have < u < /iu^. Hence, by Proposition 
14.31 and having in mind (|4.48p . wc get 

< (/ + XBfi)-^ (u) < (/ + XBp)-' ifiup) = At (/ + A^i) ' (up) <fi{I + XBp)'' (up) = ^iup. 

a 

Next we introduce the main result of this section, which paves the way for the operator Bp to generate 
an order-preserving semigroup |12j . 

Theorem 4.5. Bp is T-accretive in L^{]0,L[), and verifies the range condition 



D{Bp 



L\]0,L[)+ <zRiI + XBp) forallX>0. 



Proof. The T-accrctivity of the operator Bp is known, and that it verifies (|4.43|) also. To prove the density 

of D{Bp) in L\]0,L[)+, wc prove that VQO,L[)+ C D{Bpf ^'^°'^^\ Let < w e V{]0,L[). By K4^ . 
V G R{I + —Bp) for all ?7 e N. Thus, for each n e N, there exists m„ G D{B) such that (u„, n{v — u„)) S B. 
Since w„ = (/ + -Bp)~'^{v), by Proposition 14. 4[ we get 



(4.49) 

Let e > 0. Since 



||w„||oo <M:=M(/3,c,zy,i,||w|joo). 



n(u-'u„) = -Da(u„,0 in X>'(]0,L[), 
multiplying by w — 5e(w„), with S'^ := T^ n^y^, and integrating by parts, we get 



(v - Se{Un))n{v - Un) dx 



a.{Un,u'J{Dv - DSe{Un)) 



Then, since 



having in mind (|4.49p . we get 



-CUn{L_)S,{Un){L_) + /35,(m„)(0+). 



a(it„,u'jL)S'e(M„) > 0, 



{v - Se{Un)){v - Un) dx < 



a{un,Un)Dv 



+ -l3S,{un){0+) < -. 
n n 
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Letting e — >■ 0+, we get 

/ {v ~ UnY dx < — , 
Jo n 

and we obtain that u„ — > f in L^(]0, L[), as n — > oo. Moreover, we have u„ — > ?j in L^(]0,i[), as n -^ oo. 

Therefore v G -D(S^) and the proof of the density of DlBp) in L^{]0, L[)^ is complete. 

To finish the proof of the theorem, we only need to show that the operator Bp is closed in -L^(]0, L[) x 
L^(]0,L[). Given (w„,w„) S Bp such that w„ — > u and u„ — > f in L^(]0, L[), we need to prove that 
{u,v) e Bp. Since (u„,u„) e i5^, we have that m„ e T_Bl/+(]0, L[) and z„ := a(ii„,M^) e C([0,i]) satisfy 

(4.50) Vn^-Dzn in P'(]0,L[), 

(4.51) h{un,DT{u„)) < ZnDT{un) as measures VT G r+ 
hs{un,DT{un)) < ZnDJr'siun) as measures V 5 eV+,T e r+, 

(4.52) -z„(0) = /3 and z„(L) = -cw„(i_). 

Let T = Ta,b G 7^. Multiplying (|4.50p by T{un) and applying integration by parts (Lemma 12. ip . we get 

pL pL 

/ VnT{Un) dx = ZnDT{Un) - Z„ (L)r(u„(L_ )) - ^T(m„(0+)), 

Jo Jo 

from where it follows that 

(4.53) / z„i?r(ii„)<&(/3+||«|li)<C. 

Jo 

Here we used the boundary condition (|4.52p to be able to disregard the term related to z„(L), as it has 

the right sign. 

On the other hand, by (|4.51|) and having in mind ()2.9p . we get 

(4.54) / Zr,DT{ur,)>^ f \Di[T{un)r)\-- f T{u^fdx. 
Jo ^ Jo ^ Jo 

By ((1351) and (|I3i)) . we obtain that 

(4.55) /^ \D{[T{ur.)f)\ <- r T{u^f dx + — < ^^ + — = C. 
Jo t^ Jo c ly c 

Using the coarea formula as in the proof of Theorem 13.21 from (|4.55p we deduce that 

\DTiun)\<^ yneN. 
2a 

Then, since the total variation is semi-continuous in L^QO, L[), we have 

/ |i?r(u)|<liminf / \DT{un)\<—. 
Jo "^°° Jo 2a 

Hence, T(m) G By(]0,L[), and consequently, m e TB\/+(]0,i[). 

Since z„ = c|u„|b(u„, uJJ with |b(u„,ujj)| < 1. for all measurable subsets E c]0,L[, we have 



\Zn\ dx < C / \Un\ dx. 

IE Je 

Therefore, by Dunford-Pettis's Theorem, we can assume that 

(4.56) z„^z weakly in L1(]0,L[). 
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Moreover, since |b(M„,w^)| < 1, we also can assume that 

(4.57) b(M„, <) ^ Zf, weakly* in L°°{]0, L[). 
As ?i„ ^- u in i^(]0, L[), from (|4.56p and (|4.57p . we obtain that 

(4.58) z = cMZfe. 

As Vn -^ V in i^(]0, L[), from (|4.56p and (|4.50p . we easily deduce that 

(4.59) v = -Dz in V'{]0,L[), 
and by (l438| and (|439| . we have z € W^i^i(]0,L[) C C{[0,L]). 

Lemma 4.6. 

z(a;) = a(M(x),M'(x)) a.e.xE]0,L[ 

Proof. We use Minty-Browder's technique. Let < a < 6, let < e C^{]0,L[) and let g G C^([0,i]). 
By (|2.10p . we have that 

(4.60) / (t)[a{un,u'„)- a{u„,g')]T^i,{un){un- g)' dx>0. 

Jo 

Let us denote 

Ja{x,r) := / a{s,g'{x))ds, 



Jo 

pr pr o 

Ja'(a;,r):= / d.^[a{s,g'{x))]ds^ —{s,g\x))g"{x)ds 

and observe that 

-a(r,,b(7/„(a;)),g'(x)) [T,,,(7/„)]' = -D"^ [Ja(a:, r,,,(7/„(a;)))] + J^.{x,Ta^f,{u,,{x))), 
this we will substitute into (|4.60p . Note now that, using (|4.5ip 

4>[inD''Ta.b{lt.n) - D'' Ja{x,Ta.b{Un))] > / (j) [h{Un, DTa^bilt-n))" - D" Ja{x,Ta.b{Un))] > 0, 

/o Jo 

where the last inequality is proved using the properties of the Lagrangian (see [4]). Then we can add this 
inequality to (|4.60p : 

< / (j) [ZnDT{Un) - DJg^{x,Tafi{Un{x)))] 

Jo 

0[Ja'(a;,r(ii„(x))) - z„5'T^j,(u„)) + g'T'^^,{ur,)a{un,g')\ dx. 

'0 

Now, since 

(^L pL pLi 

(f)Zn[DTa.b{Un) - g'T^fiiUn)] = / (f) Z,iD[Ta,b{Un) - g] + (/) Z„g'(l - T^ ,,(u„)) dx 

/o 'Jo Jo 

Vn(t> {Ta,b{un) - .?) dx - / {Ta,b{un) - q) a(u„, u'J(f)' dx + / 0z„.g'(l - T^,b{un)) dx 
/o Jo Jo 

we get 

lim / (t>Zn[DTa,b{un)-g'Kt,{un)]dx<{zD{Ta,b{u)-g),4>) + \W\\oo \z\<t>{l-T'^i,{u))dx. 
"■^+°^Jo ' Jo 

On the other hand, the almost everywhere convergence of m„ implies that 

Ja' {x,Ta,b{Un{x))) ^ Ja' {x,Ta,b{u{x))) a.C. 
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and we also have (see [1], Proposition 3.13) that 

D [Ja{x,Ta^b{un{x)))] -^ D [Ja{x,Ta,b{u{x)))] wcakly as measures. 
As a consequence, we have 

hm / (j) [Ja'(a;,ra.h(u„(a;))) - DJ^{x,Ta.b{un{x))) + g'T^ ,,{un)a.{un, g')] 

n^ + oojg 'J 

i-L 



= ( Ja' {x, Ta.b{u)) - DJ^{x, T(u), (p) + 4>g' Si{u, g')T^ (,(■") d.x. 

Jo 

Consequently we obtain 

{zD{Ta,biu)-g),(b) + \\g'\\oo I \'z\4> {l ~ Ti,{u)) dx 



L 



+ / <t>a{u,g')g'T'^^{u)dx - {D[J^{x,Ta.b{u{x)))]- J^,{x,Ta.b{u{x))),^) >0 
Jo 
for all < (/) e C^QO, L[). This means that, as measures, 

Z D {Ta,b{u) -g)-D [Ja{x, Ta,b{u{x)))] + J^,{x, Ta,b{u{x))) 

+ {a(«,.9').9'T;,(^/,) + \A\WU (l - K,b{u))] C^ > 0, 
and we obtain 

z (Ta.biu) - g)' - a(M, g'KTaAu))' + s.{u,g') gX,bi^) + l^lblloo (l - T^bH) > 0. 

If a; G [a < u < b], this reduces to 

{z-a{u,g'))(u-gY>0, 

which holds for all x e ^ n [a < u < b], where C^QO, L[\Q) = 0, and all g e C^([0, L]). Being x e^n[a < 
u < b] fixed and ^ G K given, wc find g as above such that g'{x) = S,- Then 

(z(x) - a{u{x),0) (^'(a^) - C) > 0, Ve e M. 

By an application of Minty-Browdcr's method in R, these inequalities imply that 

z(x) = a(ii(a;), ^'(a;)) a.e. on [a < u < 6]. 

Since this holds for any < a < 6, we obtain the identification a.e. on the points of ]0,L[ such that 
u{x) y^ 0. Now, by our assumptions on a and (|4.58p we deduce that z(x) = a{u{x) , u' (x)) = a.e. on 
[u = 0]. The Lemma is proved. 

D 
To finish the proof we only need to show that 

-|i:)"(r(w)2)| < zD-'Tiu) as measures VT eT^, 

\D'{JseiTiu))\ < zD'JT'siu) as measures, VS* G 7^+, T e r+, 
— a(w, u')(0) = /3 and a{u,u'){L) = —cu{L^). 
These proofs are similar to those in the previous section. 

n 

From Theorem l4.5[ according to Crandall-Liggctt's Theorem (c.f., e.g., [H]), for any < uq G L^{]0, L[) 
there exists a unique mild solution u £ C([0, T]; i^(]0, L[)) of the abstract Cauchy problem 

u'{t)+Bpu{t) 3 0, u{0) = uo. 
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Moreover, u{t) = Tp{t)uo for all t > 0, where (r^(i))t>o is the semigroup m i^(]0, L[)+ generated by 
Crandall-Liggett's exponential formula, i.e., 

Tp{t)uo = lim ( /+ -Bp I uq. 

On the other hand, as the operator Bp is T-accretive we have that the comparison principle also holds for 
Tp{t), i.e., ii uo,uq £ L^{]0,L[)'^, we have the estimate 

(4.61) \\{Tp{t)uo - Ti3{t)uo)+\\i < UK - "o)+!|i. 

Obviously, by Crandall-Liggett's exponential formula, from (|4.46p . we get that for all mq G L^(]0, L[)+, 

(4.62) Tp{t){^iuo) = fiTi^{t){uo) for alH > 0. 

As a consequence of (|4.6ip and (|4.62p . for u G L°°(\0, i[) + , wc have 

0<Tp{t){u) <fj.up, with ^( = max i ^^^1^,1 I, Vf > 0. 



5. Existence and uniqueness of solutions of the parabolic problem 
This section deals with the problem 

— = {a{u,u,))^ in ]0,T[x]0,L[ 



(5.63) 



-a(M(t,0),Wa;(i,0)) ==/3 > and u(t,L) = on t e]0,T[, 
u{0,x) = uq{x) in x£]0,L[. 



To make precise our notion of solution we need to recall the following definitions given in [3]. We set 

Qt^]0,T[x]0,L[. 

It is well known (see for instance [53]) that the dual space [L^{0,T;BV{]0,L[))\ is isometric to the 
space L'^{0,T;BVQO,L[)*,BV(]0,L[)) of all weakly* measurable functions / : [0,T] -^ BV{]0,L[)*, such 
that v{f) G L°°{[0,T]), where v{f) denotes the supremum of the set {|(w,/)| : ||w||By(]o,L[) < 1} in the 
vector lattice of measurable real functions. Moreover, the duality pairing is 

{wj)^ f {w{t)J{t)) dt, 



for w e 2.1(0, T; BVQO, L[)) and / e L°°(0, T; BVQO, L[)*,BV{]0, L[)). 

By LI^{0,T,BV{]0,L[)) we denote the space of weakly measurable functions w : [0,T] — > BV{]0,L[) 
(i.e., t e [0,T] -> (w(i),(/)) is measurable for every </) G BV{]0,L[)*) such that /p \\w{t)\\ dt < oo. Observe 
that, since BVQO, L[) has a separable predual (see [1]), it follows easily that the map t G [0, T] — > ||w(i)|| is 
measurable. By Lj^^ ^{Q,T, BVQO, L[)) we denote the space of weakly measurable functions w : [0,T] ^• 
BVQO,L[) such that'the map t G [0,r] -^ \\w{t)\\ is in Ll^^QO,T[). 

Note that if w G L^{0,T; BVQO,L[)) n L°°{Qt) and z G L^{Qt) is such that there exists an element 
C G [L^{0,T;BVQO,L[))\* with D^z == ^ in V'{Qt), we define, associated with (z,^), the distribution 
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■lDxW in Qt by 

(5.64) (zDxW.ip) = —{^,ipw) — / / ■z,{t,x)w{t,x)dx^{t,x)dxdt 

Jo Jo 

for all if e V{Qt). 

Our concept of solution for the problem (|5.63l) is the following. 

Definition 5.1. A measurable function u :]0,r[x]0, L[— > R"*" is an entropy solution of ()5.63p in Qt = 
]0,T[x]0,i[ iiu & C([0,r];Li(]0,L[)), ^^O) e L^^(0,T, BF(]0,i[)) for all T e %, and z(t) := 
a(u(i),9a;M(t)) G L^{Qt), such that: 

(i) the time derivative ut of u in V^Qt) belongs to [L^(0,T; i3y(]0, L[))]* and satisfies 

pT pT pL 

(5.65) / {ut{t),ij{t))dt = - / u{t,x)e{t,x)dxdt 

Jo Jo Jo 

for all test function tp e L^(0,T; i3y(]0, L[)) compactly supported in time such that tp{t) = 

J^ eis) ds i,nd & e Ll{0,T; BV{]0, L[)) n L^iQr). 
(ii) DxZ = ut in V'{Qt), and for any w e L^{Q, T; BV(]0, L[)), the distribution zD^w defined by ([QI]) 
is a Radon measure in Qt and verifies, for all w E L^{Q,T; BVQO,L[)), the following integration 
by parts formula 

(5.66) / zDxW+{ut,w) ^ l3 w{t,0+)dt-c u{t,L_)w{t,L_) dt. 
Jqt Jo Jo 

(iii) the following inequality is satisfied 

ijhs{u,DT{u))dt+ / rihT{u,DS{u))dt< 
Qt Jqt 



JTs{u)dtri dxdt — / aiu^dxu) ■ dxij T{u)S{u) dxdt 

'Qt JQt 

for truncatures S, T G T^ and any rj G C°°{Qt) of compact support. 
In the following result we get a positive lower bound for u(t, 0+). 



Lemma 5.2. If u is an entropy solution of i5. 63\) in Qt ~ (0,r)x]0,L[, then 

(5.67) u{t, 0+) > - > 0, for almost all t e]0, T[. 

Proof. For any n G N, let Vn be the function defined by zero in ]l/n, L], 1 at a; = 0, and a straight line 
joining both values at the rest of the points. Being < </> G X'(]0,T[) fixed and taking w in (|5.66p as 
Wnit) := (t>{t)vn, we get 

(5.68) / zDwn + {ut,w„) ^13 j <f>{t) dt. 

Jqt Jo 

By ()5.65p . we have 

{ut,Wn)=- / 4)'{t) / u{t,x)Vn{x)dxdt, 

Jo Jo 

so by the Dominate Convergence Theorem, 

(5.69) lim (wt,m„) =0. 
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On the other hand, given (p £ 'D{Qt), we have 
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{zDxWn,(p) = / 0(0/ z{t,x)(p{t,x)v'„{x)dxdt. 
Jo Jo 



Hence, 
(5.70) 



z{t,x)DxW„{t,x) 



n4>{t) / z(t, x) dxdt. 
Jo 



Now, by dSUl), ((OQl) and ([5J0)) . we get 

(3 / 0(t) dt : 



hm / (f)(t)n / z(t,x) dxdt. 



Then, since |z(f,x)| < cu{t,x), by Fatou's Lemma we obtain that 



13 / 0(f) dt<c (f>{t) 



hm n / u(f, x) dx 



dt^c (t){t)u{t,0+)dt 
Jo 



from where it foUows (|5.67p . 



D 



Remark 5.3. Let u a bounded entropy solution of (|5.63p in Qt- In the proof of the next result we need 
the following time regularization. For that, given cp G 2?(]0,r[) and w G Lj^^{0,T;BV{]0,L[)), we define 
{(jywy , as the Dunford integral (see fT8\ ) 



bwY{t) ■= - 



{s)w{s) dse BV{]0,L[)**, 



that is 



((0u;)^(i),7y) = - / {cj){s)w{s),r^) ds V77 G W(]0, L[)*- 



In [2] it is shown that (0w)^ E C([0, T];BV{]0, L[)). If u is an entropy solution of (|5J3)) and p e T+, it is 
easy to see that 



\D,{M^)ym]o,L[)< 



\D.,i^{s)p{u{sm{]0,L[)ds. 



Then, by the lower-semi-continuity of the total variation respect to the L^-convergence, we have 

\D4^{t)piuitm{]0,L[) < \immi\D4Mu)rm]0,L[) 

T— >0 



<limsup- / \D,{(j>{s)p{u{s))\i]0,L[)ds. 

T^O T Jt-T 

Since the map t ^^ |I?a;((/)(i)p(u(i)))|(]0,L[) belongs to L\^^{[Q,T]), we have that almost ah t e [0,T] is a 
Lebesgue point of this map. So, for almost all t € [0, T], we have 



\D,{^is)p{u{s))m L[) ds ^' \D,m)piuitm]0, L[ 



and consequently. 



T^O 



\D.,{M^)rit)mL[) ^' \D,(mp{umi]o,L[) «-e. t. 

Respect to the existence and uniqueness of bounded entropy solutions we have the following result. 
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Theorem 5.4. For any initial datum < uq G ^°°(]0, i[) there exists a unique bounded entropy solution 
u of 15. 6S^) in Qt =]0,T[x]0,L[ for every T > such that u{0) = uq. Moreover, if u(t), u{t) are 
bounded entropy solutions of 115. 63\} in Qt =]0, r[x]0, L[ corresponding to initial data uq, uq G L°°(]0, L[)~^ , 
respectively, then 

\\{u{t)-u{t)) + \\i<\\iuo-uo)+\\i for all ^ > 0. 

In particular, we have uniqueness of bounded entropy solutions of i5.63\) . 

Proof. The comparison principle. Let 6 > a > 2e > 0, T{r) := Ta.bir) — a. We need to consider 
truncature functions of the form Sed{r) := T^{r - l)+ = Tij+e{r) - I G T+, and Sl{r) := T^{r - l)~ + e ^ 
Ti_^^i{r) +e-l€ T+ , where I > 0. Observe that Sl{r) = -T,{1 - r)+ + e. Let us denote 

4eii^)^ [T{s)T,{s-l)+ds, 
Jo 

JTei{^)= T{s)T,{s-l)-ds = - I T{s)TS - s)+ ds. 



JQ JO 

Then, JTS,.iir) = J^^^jir) and Jrsiir) = -^r,e,;('') +e-'^r(r). 

Let u, Ti be two entropy solutions of (j5.63p corresponding to the initial conditions uq, uq G (^^(]0, i[)) , 
respectively. Then, if z(i) := a{u{t) , dxu{t)) , z{t) := a{u{t) , dxu{t)) , and li,l2 > e, we have 

/ J+ ^ i^iu{t))dtvit) dxdt 
Jo 

(5.71) +/ / vit)[hTiuit),DxS,,i,{u{t))) + hs,,M*),DxTiuit)))]dt 

Jo Jo 

/ z(t)9,77(t) T{u{t))S,,i,{u{t)) dxdt < 0, 

"'0 

and 

/ >/j^e ia''"^^))^*'^ "^^"^^ ^ ^ / / JT{u{i))dtrj{t) dxdt 
Jo ' ' ^ Jo Jo 

(5.72) +f f ,j{t)[hT{u{t),DxSHu{t))) + h^i,{u{t),DxT{u{tmdt 

Jo Jo 

/ z{t)dxr]{t) T{u{t))Si^iu{t)) dxdt < 0, 
Jo 

for all ?7 e C°°{Qt), with 77 > 0, r]{t, x) = (j){t)p{x), being e X>(]0, T[), p e P^O, L[). 

We choose two different pairs of variables (t, x), {s,y) and consider u, z as functions in {t,x), u, z in 
{s,y). Let < ^ S I?(]0,T[), f/; g I?(]0, L[), p^ and p„ sequences of moUificr in M. Define 

, , , ._, . / 1 + s\ / X + y 

r]m.nV;X,s,y) := Pm[x - y)pn(t- s)(p —-— V 
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For (s, y) fixed, if we take in ()5.7ip h = u(s, y), we get 



+ 



"'0 
T ^L 



(5-73) +/ ^ 77,„,„[/iT(w(t, a;),£'a;5'e,ji(s,y)(u(i,a;))) + /is^_(^_^j(M(i, x),I?a;r(u(i,.T)))] di 

z{t, x)dxr]m,n T{u{t,x)) S^ji(s,y){.u{t,x)) dxdt < 0. 



^0 
T ,.L 



Jo 



Similarly, for (i,x) fixed, if we take in (j5.72p I2 = u{t,x) we get 



"^Te«(ta;)("(*'y))^s??m,n rfyds- e / / JT{u{s,y))dsrim,n dyds 

a Jo ' Jo Jo 

(5-74) +/ / VrnAhT{^s,y),DyS:^^'^''\u{s,y))) + hg.,,..,{u{s,y),DyT{u{s,ymds 

Jo Jo 

T nL 

z(s,y)a.y77„,,„ T{u{s,y)) S:^^'-'^\u{s,y)) dyds < 0. 
"'0 



Wc integrate (|5.73p in {s,y), (|5.74p in {t,x), and add the two inequalities. Using that a > 2e, and since 



r]m,nhs,,^(,_y-,iu{t,x),DxT{u{t,x))) dsdtdy > 

QtxQt 



and 



r]rn,nhg^(t.^){u(s,y),DyT{u{s,y))) dsdtdx > 0, 

QtxQt 
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we get 






T,e,«(s,y)('"(^'^))^*'?™." + '^T,e,«(t,2^)('^(*'2^))^s^™>» ) dsdtdydx 



-el JT{u{s,y))dsr]m,n dsdtdydx 

JtxQt 

'nm,.nhT{u{t,x),DxS^,u{s,y){u{t,x))) dsdtdy 

■tX-Qt 

Vni.nhT^uis, y), DyS'^'^*''-''^ (u(s, y))) dsdtdx 
Qry-Qr 

^^■^^^ -/ z{s,y)d^Ji^.nT{u{s,y))Sf'^-\u{s,y)) dsdtdydx 

z{t, x)dyr]jn,nT{u(t, x))S^^u{s,y)iu(t, x)) dsdtdydx 

QtxQt 

T+iuit, x) - u{s, y))[T{u{t, x))z{t, x) - T{uis, y))z(s, y)] 

QtxQt 

xidxVm.n + dyrim,n) dsdtdydx 
-e / T{u{s, y))z{s, y)(9a;?7,„,„ + dyijm.n) dsdtdydx < 0. 

'QtxQt 

Let I2 be the sum of the third up to the sixth terms of the above inequahty. From now on, since u, z are 
always functions of (i, x), and u, z are always functions of (s, y), to make our expression shorter, we shall 
omit the arguments except when they appear as sub-index and in some additional cases where we find it 
useful to remind them. We also omit the differentials of the integrals. 

Working as in the proof of uniqueness of Theorem 3 in g], we obtain that i/2 > ll</'lloollV'l|ooo(e). Hence, 
by ([5T75)) . it follows that 



txQt ^ ' ' ' 

T^{u - u)[T{u)z - T{u)z\{dx'nm,n + dy-qm,n) 



Qt-xQi 

-e / T{u)z{dxilm,n + 9y7?m,„) < eo(e) + e / JT{u)ds'nm,n- 

^txQt JQtXQt 

Then, dividing by e and letting e — )• we get 

QtxQt ^ ^ 

sign|]'(M - u)[T{u)z - T{u)z]{dx11m,n + dyrim^n) 
QtxQt 

T{u)z{dx-qm.,n+dy'qrn,n) < / JT{u)ds'ri,n,n 

QtxQt JQtxQt 
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where 

4-9"/0 = r T{r')sign+{r' - l)dr' / e M, r > 
Jo 
and 

JT,s^gn^r) = j T{r')sign^ {r' - l)dr' Z G M, r > 0. 
Jo 
Now, letting to — >■ oo, we obtain 

rT .T .L 

/ / sign^ {u{t, x) -u{s,x))[T{u{t, x))z{t,x) -T{u{s,x))z{s,x)]dxXn 
'0 Jo Jo 

i^T t^T pL t^T pT pL 

/ / T{u{s,x))z{s,x)dxXn < / / / JT{u{s,x))dsXn 
/o Jo Jq Jo Jo Jo 

where X„(i, s,x) := pn(t — s)(J){^^)iIj{x). We set ^ — tlj^. ^ 2?(]0, L[) t ^]o.l[ in the last expression and 

taking limit as fc — > +00, we have 

rT rT rL _ _ 

■'0 ■'0 

I'T pi pL 

+ lim / / / Kn{t,s)sig-n:^{u{t,x) ~u{s,x))T{u{t,x))z{t,x)dx'>pk{x) 

k^+oojg Jo Jo 
pT pi pL 

(5.76) - lim / / / K„(i, s)signj(w(i, x) -u(s, a;))r(u(s, a;))z(s,x))52;'(/'fc(a;) 

fc^+ooJo Jo Jo 

pi pi pL 

+ lim / / / K„{t,s)T{u{s,x)z{s,x))dx^k{x) 

fc^ + ooJo Jo Jo 

pi pT pL 
-ill JT{u{s,x))dsK.nit,s), 

Jo Jo Jo 
where K„(i, s) := /5„(i - s)0(^). 

Let us study the second and the third term of the above expression. Let 

pi pi pL 

Ik ■= / / Kn{t, s)sign^ {u{t, x) -u{s,x))T{u{t,x))z{t,x)dxi'k{x) 

Jo Jo Jo 

T i-T pL 

/ / K„{t,s)sign^{u{t,x)-u{s,x))T{u{t,x))z{t,x)dx{ipk{x)-l). 

/O "'0 Jo 

Let Hn{s,r) := K„(r, s)signo (u(r) — u{s))T{u{r)). For r > 0, we define the function (k„(s))'^, as the 
Dunford integral (see Remark [53| 

{K„{s)y{t):^-J H,,{s,r)dr. 
Then, 

/•T^ z*^ r^ 
Ik = lim / / (K„(s))^(t)z(t,a;)92;[^fc(a;) - l]rfa;(itds 
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pT pi pL 

= -lim/ / / [i^k{x)~lHt,x)D,,{{K,,{s)y{t)))dsdt 



lim / (m(, (K„(s))'^(?/>fc(x) - 1)) ds 



+clim / / u{t,L\{K.Js)Y{t){L^)dtds 
^^oJo Jo 

-P lim / / (^„(s))-(i)(0+) dtds = II + II + II + It. 
Notice that 

ll^cl / u(i,i_)K„(i,s)sigii([(w(t,L_) -u(s,L_))T(M(i,i_))dfds 
Jo Jo 

and 

/4 = -/?/" f K,,{t,s)sign+{u{t,0+)-u{s,0+))T{u{t,0+))dtds. 
Jo Jo 

By Remark l5.3i we get 

(5.77) |Z?,(K(s))^(t))|(]0,L[) ^-^ \D.,iK„it,s)sign+{u{t)-u{s))T{u{tm{]0,L[). 



Using (|5.77p . we get 



pT pT pL 

\ll\<c\\u\\L^iQT) / / {l-Mx))\DA^n{t,s)sign+{u{t)~u{s))T(u(tmdtds, 

Jo Jo Jo 

which implies limfc_^oo ll = 0. Let us deal with I^. We have 

II = lim r r r u{t,x) ""^''^^^^ ' ""^''^\ Mx) - l)dxdtds. 
■r^oJo Jo Jo T 

Let 

qir) := sign+(T - u{s, x))T{t), Q{r) := / ^(t) dr. 

Jo 

Since 5 is non-decreasing, (3(r)—(3(r) < g(r)(r—r). Then, changing variables, since -ff,i(s,t) ~ q{u{t))Kn{t, s) 

"-^"Jo Jo Jo T 

fT fT rL ^Q{u{t,x))-Q(u(t-T,x)) , , , 

> lim / / / (l--0fc(a:))K„ft, s) ^^ ^ ' ^^ ^^^^ '-^^ dxdtds 

^^oJo Jo Jo T 

(5.78) 

= lim r r [\l -~ Mx))Q{u{t, x)) ''"^^' '^ - ^"(^ + r, s) ^^^^^^ 

^^OJo Jo Jo T 

/ / {I- i^k{x))Q{u{t,x))dtKn{t,s)dxdtds, 
Jo Jo 
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from where it follows that limfc^oo I^ > 0- Taking into account the above facts, we get 

lim Ifc > -/3 / / Knit, s)sign^ {u{t, 0+) - u{s, 0+))T{u{t, 0+)) dtds 
(5.79) 

+c / / u(t, L^)Knit, s)aignQ {u{t,L^) — u{s,L^j)T{u{t,L^j) dtds. 
Jo Jo 

Working similarly, we obtain 

pi PI pL 

lim / / / Kn{t, s)sign^ {u{t, x) ~u{s,x))T{u{s,x))z{s,x)d:j:i:k{x) 



k—^oc 



Jo Jo 



(5.80) > /3 /" f Kn{t, s)sign+(u(i, 0+) - u(s, 0+))T(u(s, 0+)) dtds 

Jo Jo 

— c / / u{s, L ^) Kn{t, s)signQ {u{t,L^) — u{s, L^))T{u{s,L^)) dtds. 
Jo Jo 

Analogously, 

pi pi pL 

lim / / / Kn{t,s)T{u{s,x))z{s,x)dxipk{x) 
(5.81) 

f^ r - - r r - 

>c / u{s,L-)Kn{t,s)T{u{s,L-)) dtds — /3 / / Kn{t,s)T{u{s,0+))dtds. 

Jo Jo Jo Jo 

From ((5J6)) . by (i5J9| . ([STSO]) and dSH]), we have 

J J (^T,sig„,«(.,:.)("(*'^))^*'*"(^''S) + '^T,sign,n(t,:.)("(«'2;))^«'^»(*'«)) ^tdsdx 

+c / u{t,L^)Kn{t,s)signQ{u{t,L^) — u{s,L-))T{u{t,L^))dtds 

Jo Jo 

— c / / u{s, L -) Kn{t, s)signQ {u{t,L^) — u{s,L^))T{u{s, L^)) dtds 
Jo Jo 

(5.82) -/3 f f Knit, s)sign+{uit,0+)^u{s,0+))T{u{t,0+)) dtds 

Jo Jo 

+/3 / / Kn{t,s)sign+{u{t,0+)~u{s,0+))T(u{s,0+))dtds 
Jo Jo 

r r- - r r - 

+c / u{s,L-)Kn{t,s)T{u{s,L-)) dtds — j3 / Kn{t,s)T{u{s,0+)) dtds 

Jo Jo Jo Jo 

pi pi pL 

< JT{u{s,x))dsKn{t,s) dtdsdx. 

Jo Jo Jo 
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By Lemma 15.21 we have 

(5.83) u(t, 0+) > - > 0, u{s, 0+) > - > for almost everyt, s > 0. 

c c 

Letting a ^- 0, dividing by b and letting & ^- in (|5.82|) . we obtain, 

rT rT pL 

{u{t, x) — u(s, x))'^ {dtHn{t^ s) + dsKn{t, s)) dtdsdx 
10 Jo Jo 

+c / / u{t, L _) Kn{t, s)sign^ {u{t, L _) — u{s,L-))signQ {u{t,L^)) dtds 
Jo Jo 

—c / u{s, L -) Kn{t, s)sign.Q {u{t,L-) — u{s,L-))sign.Q {u{s,L-)) dtds 

Jo Jo 

~P / K„(t, s)signj(u(t, 0+) -M(s,0+))[signo(u(i,0+) -signo(u(s, 0+))]dids 

Jo Jo 

+c / / u{s, L-)Kn{t, s)signQ {u{s,L^)) dtds ~ /3 / / Kn{t,s) dtds 
Jo Jo Jo Jo 

t^i pi pL 

< / / u{s,x)dsKn{t,s) dtdsdx. 

Jo Jo Jo 
Having in mind (|5.83p . the fourth term of the above expression vanishes. Moreover, the sum of the second 
and third term is non-negative. On the other hand, since Us = Dx (z) in the sense given in (ii) of Definition 

EH 

i-T pT pL pT 

u{s,x)dsKn{t,s)dxdtds = — / {Us,Kn{-,t)) dt 
/o Jo Jo Jo 

i-T rT j-T j-T 

u{s, L-)Kn{t,s) dtds — /3 / / K„{t,s) dtds. 
Jo Jo Jo Jo 

Therefore, 

rT rT rL 

{u{t, x) — u{s, x))'^ {dtKn{t, s) + dsKn{t, s)) dtdsdx < 0. 
/o Jo Jo 
Letting n — > cx), 

T rL 



{u{t, x) - u{t, a;)) + 0'(t) dxdt < 0. 
/o Jo 
Since this is true for all < £ 2?(]0, T[), we have 

d f^ 

— / (u(t, x) - u{t, x))+ dx < 0. 
dt Jo 

Hence 

rL pL 

{u{t,x) — u{t,x))^ dx < / {uq{x) — uo{x))^ dx for all < > 0, 
'0 Jo 

which finishes the uniqueness part. 

Existence of bounded entropy solution. Given < Uq G L^{]0,L[), let u{t) ~ Tp{t)uQ, being (r^(t))t>o the 
semigroup in L^(]0,L[)~'" generated by the accretive operator Bp. Then, according to the general theory 
of nonlinear semigroups (|12j). we have that u(t) is a mild-solution of the abstract Cauchy problem 

u'it)+Bpu{t)^0, u{0)^uq. 

Let us prove that, assuming < wo S L°°{]0, L[), then w is a bounded entropy solution of (|5.63p in Qt. 
We divide the proof of existence in several steps. 
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Step 1. Approximation with Crandall-Ligget's scheme. Let T > 0, A' > 1, At = -^, tn ~ nAt, n = 
0, . . . ,K. We define inductively u""*"^, n = 0, . . . , A' — 1 to be the unique entropy solution of 



m"+i - u''' 



(5.84) 

where u^ = 7^0 • 
If we define 



(a(w"+\(u"+i)'))' = in ]0,L[ 



At 
-a(u"+i(0),(u"+i)'(0)) =^>0 and w"+i(A_) = 0, 



K-l 



I'^it) := u%o,t,m + E ""^l*.,*.+i]W' 



by Crandall-Liggett's Theorem, we get that u^ converges uniformly to u e C([0, T], i^(]0, A[)), as K — )- oo. 
We also define 

^-1 ,,n+l _ „,n 
n=0 

and 

K-l 

z^(t) := a(^\ (^i)')X[o,t,](t) + E ^K^'' K+')')X]*„,Wi] W- 

Since w"+^ is the entropy solution of (|5.84p . we have 

(5.85) ^^ (t) = D^z^ (t) in V'{]0,L[), Vte]0,r] 

(5.86) z^{t){L) = -cu^{t + At){L^), Vt e]0,T - At], -z^^'(t)(0) = /3, Vt <E [0,T] 
and for all S* G 7^+, T e r+, we have Vt g]0, T - At] 

(5.87) h{u^{t + At),D^T{u^{t + At)) < z'^{t)D.^T{u^{t + At)) as measures 

hsiu^it + At), D^T(u^(t + At)) < z^{t)D^JT's{u^{t + At)) as measures. 
Note that (|5.87p is equivalent to 

||D^((T(u^(t + At)))2)| < z^{t)DlT{u^{t + At)) as measures. 

Since a(w"+i, {u"+^)')D^T{u"+^) > h{u'^+^,D^T{u''+^)) as measures in ]0, L[, using (gH) we can write 

/i(u"+\7^,r(7."+i)) = a(«"+i,(u"+i)')(T(«"+i))'£i + ^|i?^[(T(u"+i))2]| 

> ^\iiT{u-+'))'y\C' - ''^{T{u-^^)fC} + \\D%{{T(u-^^)f\\ 

= 'L\DA{T{u-^^)f\\--{T{u-^^)fL\ 
2 J/ 

from where we get the following inequality as measures 

(5.88) z'<{t)D,T{u''it + At)) > ^|i?,[(T(u^(t + At)))2]| _ -(T{u'<it + At)))\ 

Lemma 5.5. There exists M := AI{P,c,i>, L, ||uo||oo) such that 

(5.89) \\u^it)\\oo < M VA' e N and Vt e [0,T]. 
Consequently, \\u{t)\\oo < M Vt e [0,r]. 
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Proof. Since 



ioi n ^ 0, . . . ,K — 1, 



by Proposition 14.41 if /i := niax{ ^^^"g^^°° , 1}, we have 

0<u^ = {I + AtBp)~'^ (wo) < fJ-up. 
Then, repeating this process, we obtain 

< u"+^ = {I + AtBpy^ (u") < (/ + AtS^)"^ {fiufs) 

= n(^I + AtBji^ (up) <n{I + AtBp)~^ (w^g) = Mu^, 
and the proof concludes. 



D 



Step 2. By (jSlSQ)) . ||z^^(t)||oo < C for all T^f e N and a.e. t G [0,T]. Then we may assume that z^^' ^ 
zgL'^{Qt) weakly*. Moreover, since z^'(t) = cu'^(i + Ai)b(M^(t + Ai),5^M^(t + At)) V<e]0,r-Ai], 
with \\h{u^{t + At),dxU^{t + Ai))||oo < 1 and u^ converges uniformly to u in C{[0,T],L'^(]0,L[)), we 
may also assume that h{u^ {t + At), dxU^ {t + At)) -^ Zb(t) G L°°{Qt) weakly* and 

(5.90) z{t) = cu(t)zb(t) for almost ah t e [0,T]. 

Given w G BV{]0,L[), from (ISTSSj) and I^Mi), it follows that for each t e]0,T] 



S, {t,x)w{x)dx 



z"^ {t)Dw + z"^ {L)w{L) + /3^(;(0^ 



< C\\w\\bv(]O.L[) + |z^'^(i)u'(i)l <{C + Cfl\\up\\oo)\\w\\BV{]O.L[), 

where the continuous injection oi BVQO, L[) into L°°{]0, L[) was used. Thus, ||C^(i)||_BV'(]o.L[)' < C, V K E 
N and t €]0,r]. Consequently, {^^} is a bounded sequence in L'^{0,T; BV{]0,L[)*). Now, since 
L°°(0,T;B1/(]0,L[)*) is a vector subspace of the dual space {L^{0,T; BVQO,L[)))\ we can find a subnet 
^" of ^^ such that 

r^ee (i'(0,r;W(]0,L[)))* weakly*. 

Working as in the proof of Theorem 5.5 of |4], Step 2, we can prove that (|5.65[) holds and ut ~ 
D.,z mV'iQr). 

Step 3. Next, we prove that ut = D^ z in the sense given in (ii) of Definition 15.11 To do this, let us first 
observe that we can prove, as in the proof of Theorem 5.5 of [4], Step 4, that the distribution zDw in Qt 
defined by (|5.64p is a Radon measure in Qt for all w G L^(0, T; BV{\Q, i[)), and also that 

{zDxW, ip) =lini / z"[t,x)Dx'w[t,x)Lpx[t,x)dxdt. 

" Jo Jo 

From where it follows, combining with (|5.86p and integrating by parts, 

pi pL PI pL 

zDxW = lim / / z"{t)Dxw{t) dt = - lim / / w{t, x)DxZ°'{t, x) dxdt 
"^ Jo Jo " Jo Jo 



Qt 



-lim 



z"'{t,L)w{t,L^)dt- / z°'{0)w{t,0+)dt 
Jo 



lim 



-(C",u')-c/ u°'{t + At){L^)w{t,L^)dt + p w{t,0+)dt 
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= -{ut,w)-c u{t){L^)'w{t,L.)dt + l3 w{t,0+)dt, 
Jo Jo 

and (|5.66p holds. 

Step 4. Let T = Ta,b be any cut-off function, let j be the primitive of T. Let < ^ £ r»(]0, T[) . Multiplying 
(|5.84p by r(u"+^)(/)(t), t G (i„, i„+i] integrating in (t„, i„+i] x]0, L[ and adding from n = to n = K— 1, 
we have 

^~l rtn+i pL „.n+l _ n i-T pL 

Y. m XT T{u^+^)dxdt+ m z''{t)D.,{T{u''{t + At)))dt 

(5.91) 

= J2 / (/3</'(t)TK+i(0+)) - c</.(i)M"+i(L_)rK+i(L_))) dt. 

Since <?!) has compact support in time in (0,r), for K large enough, performing like in (|5.78p . we have 



n=0 



Hence, from (|5.9ip it follows that 

z^{t)(t){t)D.^T{u^'{t + At)) dt 
la JQ 
(5.92) 

- f f ^•("''W) '^^^^'At~^^^ dxdt + f p^{t)T{u^{t + At, 0+)) dt. 

Given e > 0, if we take into ([02]) any test < e ^(JO, T[) such that 0(t) = 1 for f e]e, T - e[, having in 
mind (|5.89|) . we get 

rT-e pL 



/ z^{t)D.^T{u^{t + At))dt 



< 




T pL 



Jiu^'it)) '^^^^ y ^^^ da-dt+ / /?T(7.^(t + At,0+))dt<C. 



"'0 ^^ "'0 



This implies that {z^{t)Dx{T{u^{t + At)))} is a bounded sequence in 

^\oc Mi(0' 2^1 -^(]0i -^[))> where A^(]0, i[) denote the space of bounded Radon measures in ]0, L[. 



On the other hand, by ([57 



i-T-e i-L 



/ / z^{t)D^T{u^{t + At))dt 



-I J V \D.[{T{u"{t + At)))^]\dt~ j Y ^(T(^^^(t + At)))2dt. 
Hence 

/ |i?,[(r(u^(t + At)))2]|dt<f^ + f2il£^ = C, 
where by the coarca formula it follows that 

/ |£)^r(w^(t + At))|dt <C. 



\D,T{u{t))\ < liminf / \D,T{u'' {t))l 

K—^oc 
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Moreover, by Lemma 5 of [2], the map t H' \\T{u^{t))\\Bv{]o,L[) is measurable, then by Fatou's Lemma 
and ()5.93p . it fohows that 

pi — e pL pi — € pL 

(5.94) / hminf / \D ^T [u^ {t + M))\ dt < VmiiTLi / \D ^T {u^ [t + M))\ dt < C . 
Now, since the total variation is lower semi-continuous in L^(\0,L[), we have 

/o ■^^'^ ^0 

thus we deduce that T{u{t)) e BVi]Q,L[) for almost all t G (0,r) and consequently u{t) G TBV+{]Q,L[}. 
Then, by (|5.94p . applying again Lemma 5 of [5], we obtain that 

(5.95) T{u{-))eLl,^^{Q,T,BV{]0,L[)). 

Step 5 (identification of the field). Let us now prove that 

(5.96) z(t) =a(M(t),5^u(i)) a.e.ie]0,r[. 

Let < e V{Qt) and g e C'^{[0,L]). Assume that (j) ^ ■q{t)p{x) with -q G V{]{),T[) and p e V{]0,L[). 
Let < a < fe, and T = Ta.b- Let j denote the primitive of T. Recall that 

Ja{x,r)= / a.{s,g'{x))ds and Ja'(a;,r)= / dx[a.{s , g' {x))] ds 
Jo Jo 

For simplicity, we write 

D2MX, T(u^(t + At))) := D, [J^{x, T(w^(t + Ai)))] - J^,{x, T{u'^{t + At))). 
Working as in the proof of step 6 of Theorem 3 in [4] we find out that 

(5.97) [i?2 Ja(x,r(u^(i + At)))r' ^ a.iu'^it + At),g')d,[T{u''{t + At))]. 
Using ([CT7)) . (PJU|) and ([ESS)) , we obtain 

T rL 



'0 JO 



(j) z^^(t)AT {T{u^{t + At))- g) dt 

I I cf,[D2Mx,T{u''{t + At)))-ei{u''{t + At),g')g']dt 
Jo Jo 

(j) [■z^{t)D^Tiu'^it + At)) - z^{t)g' + &{u'^ {t + At),g')g'] dt 



Jo 

T rL 



Jo 

I 4>{[D2Ux,T{u''{t + At)))r + [D2Je.{x,T{u''{t + At)))Y] dt 

JO 

T pL 

(f> {a{u'^{t + At),g') - z^(t)) {g' - d^T{u^{t + At))) dt 

(p [i.'^{t)D%T{u^{t + At)) - [D2J^{x,T{u'^{t + At)))]''^] dt 
Jo 

> f f ^ [z''{t)D'^T{u''{t + At))) - [D2Mx,T{u'^{t + At)))]'] dt 
Jo Jo 

>ff cj,\^\D'jT{u''it + At))')\-[D2Ux,Tiu''{t + At)))]' dt. 
Jo Jo '-^ 
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Again working as in the proof of step 6 of Theorem 3 in [4] , we get 



T rL 



'0 Jo 
Therefore, we obtain 

i-T rL 



- \D^,iT{u^it + At))2)| - [D,i.Ux, T{u^{t + At)))]' 



dt> 0. 



(5.98) 



,K 



(j) z^ {t)D^{T{u'^ {t + At)) - g) dt 



Jo 

T i-L 



f f (b[D2Ux,Tiu''{t + At)))-eL{u''{t + At),g')g']dt>0. 
Jo Jo 



Now we shah bound from above the first term. By (j5.85p and for At small enough, performing like in 
(I^TTg]) . we get 

rT rL r-T pL 

(P{t,x)T{u^{t + At))D.^z^\t)dt^ / 4){t,x)T{u^{t + At)^^{t)dxdt 



JO 

T rL 



JO 



> 



Jo 



t^L^^h^^Z^^^^j^u'<it))dtdx 



Then, integrating by parts, we have 

rT rL 



(j){t) z" {t)D^{T{u'^ {t + At)) - g)dt<- 
Jo Jo Jo 



^ '■^0(i-Ai)-0(t)., K 



At 



j{u'^{t))dtdx 



T rL rT i-L 

<l>{t)g£,'^{t)dtdx- / da:<l>{t) z'^{t)[T{u'^{t + At))~g]dtdx. 

Jo Jo JO 



Thanks to this inequality we arrive from (j5.98p to 

_ r 1^ Ht - At) - 0(t) 

^0 ^0 



At 



j{u^(t))dtdx^ 



T rL 



0(05? (t) dtdx 



(5.99) 



T rL 



JO 
T rL 



Jo 



d^(p{t) z^(t)[r(u^(t + At)) - 5] dtdx 

4>{t) [D2J^{x,T{u^{t + At))) - SL{u^{t + At),.g')5'] dt > 0. 



Letting /-iT — > cx) in (|5.99p and having in mind that 

-D2 Ja(a;, T{u'^{t + At))) -^ D2Ja{x, T{u{t))) weakly as measures 
we obtain 

pi pL pi pL 

/ / dt(l){t)j{u{t))dt + {ut,H) - / [T{u{t))~g]z{t)d:,(j){t)dxdt 

Jo Jo Jo Jo 

(5.100) 

+ / / <P{t) [-D2Mx,T{u{t))) + a{u{t),g')g'] dt > 0. 
JO Jo 

By (EH, 



pT pL pT pL 

{ut,<t>g)=-l I z{t)gdx(f>{t) dtdx - / z{t)g'(f>{t) dtdx 

Jo Jo JO JO 
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and we can rearrange ()5.100p in the following way 

dt(f>{t)j{u{t))dtdx- / z{t)g'4>{t)dtdx- / T{u{t))z{t)dx(j){t) dxdt 

In JQ Jo Jo Jo Jo 

(5.101) 

4>{t) [-D2Mx,T{u{t)))+a{u{t),g')g'] dt > 0. 
/o Jo 

Now, for T small enough and using again the trick in (|5.78p . we have 

dt(j){t,x)j{u{t,x))dxdt= Mm / ^^^^^—-^—^^j{u{t,x))p{x)dxdt 

Jo '"^OJo Jo -T 

< lim / / u(t,x)p(x)—-(r/T(u)Y(t,x)dxdt, 
'"^oJo Jo dt 

where we used again the notion of Dunford integral (see Reniark l5.3p . Using (|5.65p . we have 

u{t)p-hT{u)nt)dxdt = -{ut,p{TjT{u)y{-)) 
Jo "-^ 

= - \im{C,p{vT{u)y{-)) = - lim / (Z?,z"(i), p- / v{s)T{u{s)) ds) dt 

"'Jo T Jt-T 

lim/ / z"'{t)D.Jp- I ri{s)T{u{s))ds] dt^ lim I [ d^-pz"{t)l -t]{s)T{u{s)) dsdxdt 
"^ Jo Jo V T Jt^T J '^ Jo Jo Jt-T T 

+ lim / / pz°'{t)D^[{r]T{u)y{t)]dt= - rjis) T{u{s))z{t)d.^pdxdsdt 

" Jo Jo Jo '^ Jt-T Jo 

pi pL pi pL 

+ lim / / pz''{t)d4{ijT{u)y{t)]dxdt + \im / p z'' (t) D'J{7]T {u)y (t)] dt 

" Jo Jo " Jo Jo 

< - V{s) / T{u{s))z{t)d^pdxdsdt + - vis) / pz{t)d^[T{u{s)))dxdsdt 

Jo '^ Jt-T Jo Jo '^ Jt-T Jo 

T ^ rt rL 



^ '•^ d 



T]{s) / cMp\D^^[T{u{s)]\dsdt. 
''' Jt-T Jo 



Taking limits when r — > 0, having in mind (j5.89p . we obtain 



T pL pT pL 

/ dt(j>it)j{u{t))dxdt< ri{t) T{u{t))z{t)d^pdxdt 
Jo Jo Jo 

pi pL pT pL 

+ T]{t) pz{t)d.^T{u{t))dxdt + cM v{t) p\D'jT{u{t)]\dt. 
Jo Jo Jo Jo 

From (|5.10ip . all gathered together reads 

pi pL pi pL pi pL 

0<-/ / (l){t)z{t)g' dxdt + r]{t) p z{t)d^{T{u{t))) dxdt + cM rj{t) I p\DlT [u{t))\ dt 
Jo Jo Jo Jo Jo Jo 

T pL 

/ 4> [-D2MX, T{u{t))) + aL{u{t),g')g'] dt. 
Jo 



ON A NONLINEAR FLUX-LIMITED EQUATION ARISING IN THE TRANSPORT OF MORPHOGENS 39 

Using that D2Mx,T{u{t)) ^ aL{u{t),g')dx{T{u{t))) + [D2Mx,T{u{t)))Y , this is written as 

pT pL pT pL 

0<cM ri{t) p\DlT{u{t))\dt- / (t>[D2J^{x,T{u{t))Y dt 

Jo Jo Jo Jo 

+ f I y- d.{T{u{t)))][a.{u{t),g') - z{t)]<l>dxdt. 
Jo Jo 
As measures, 

cM\D',T{um - [D2Ux,T{u{t))r + [g' - d.,{T{u{t)))Mu{t),g') - z{t)]C'- > 0. 

Taking the absolutely continuous part and particularizing to points .t G [a < u{t) < b], this reduces to 

[g'~dMt)Muit),9')-<t)]>Q, 
an inequality which holds for all (i, a;) € Sr\[a <u< 6], where S* C]0,T[x]0,L[is such that /;2(]0,r[x]0,L[\S') = 
0, and all g G C^([0, L]). Being {t, x) e S Ci [a < u < b] fixed and ^ S M given, we can find a function g as 
above such that g'{x) = ^. Then 

(z(<, x) - a(M(t), 0) {d^u{t, x)-^)>0, VC e K and V(t, x) e S r][a < u < b]. 

By an application of Minty-Browdcr's method in M, these inequalities imply that 

z{x) = a(u(t, x), dxu{t, x)) a.e. otiQt n [a < u < 6]. 

Since this holds for any < a < 5, we obtain (j5.96|) a.e. on the points of Qt such that u{t,x) ^ 0. Now, 
by our assumptions on a and (|5.90|) we deduce that z[x) = a(u(a;), u'(a;)) = a.e. on [u = 0]. We have 
proved ((05|) . 

Step 6. The entropy inequality. Given 5* £ V^ , T G T^ and (f> e VIQt), working as in the proof of (|5.92p 
wc can get 

rT pL 

<f> z'^{t)D.x (r(u^(f + At))S{u^{t + At))) dt 

(5.102) <J^ J^ Jrs{u-{t)) ^^'^-f^-^'^ dxdt 

/ z^{t)dx(t)T{u^{t + M))S{u''^{t + At)) dxdt 
Jo 

and the fact that {z'^{t)Dx (r(u^(t + At))S'(w^(t + At)))} is abounded sequence in ii„^(0,T;X(]0, L[)). 
From here, as in the proof of Theorem 4.5 of [6], we can get that the sequences {z^ [t)DxJT' s[u^ [i + ^i))} 
and {z^{t)DxJs'T{u^{t + At))} are bounded in L]^^{Q,T;M{]Q,L[)). This allows us to define, up to 
subsequence, the objects /i|^, ^^ S M[Qt) by means of 

((/), A^i) = lim / j d^z^{t)DxJT's{u'^{t + M))dt, yc^eCdQT), 
K Jo Jo 

(cj), /x|) = lim / / z''{t)D,Js'T{u''{t + At)) dt, V0 e C^Qt)- 
^ Jo Jo 
Then, passing to the limit in (|5.102p . wc obtain 

((/),/i^) + (0,M|) < / / JTs{u{t))dt(l){t)dxdt 
Jo Jo 
(5.103) 

/•T pL 

z{t)dx(t)T{u{t))S{u{t)) dxdt, Vcf) e V{Qt), 



Jo 

T i-L 
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Working as in proof of Lemma 4.11 in [6], we can get the following result. 
Lemma 5.6. For S G r+,T e T+ , we have that n^ > hs{u,DT{u)). 

By the above lemma and (|5.103p we obtain the entropy inequality 

/ / (l)hs{u,DT{u))dt+ / (j)hT{u,DS{u))dt< 

Jo Jo Jo Jo 

/ / Jts{'u-)4'' dxdt — / / a.{u,dxu) ■ dx4) T{u)S{u) dxdt 
Jo Jo Jo Jo 

for truncatures S G V^ , T € T^ and any smooth function (p of compact support. 
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